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Determination of the anomalous dimension of gluonic operators 



Abstract. Using large Nf methods we compute the anomalous dimension of the predominantly 
gluonic flavour singlet twist-2 composite operator which arises in the operator product expansion 
used in deep inelastic scattering. We obtain a d-dimensional expression for it which depends on 
the operator moment n. Its expansion in powers of e = (4 — d) /2 agrees with the explicit exact 
three loop MS results available for n < 8 and allows us to determine some new information on the 
explicit n-dependence of the three and higher order coefficients. In particular the n-dependence 
of the three loop anomalous dimension 7 S9 (a) is determined in the 02(G) sector at 0(1/Nf). 
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1 Introduction. 



Perturbative quantum chromodynamics, (QCD), is widely accepted as being the tool to describe 
high energy scattering processes involving the strong interactions. (For a review see, for example, 
]l]].) This is primarily due to its property of asymptotic freedom which implies that as higher 
energies are attained the strength of the strong coupling decreases and perturbation theory 
becomes a more valid approximation. One of the current topics of interest is in constructing 
the perturbative series at a higher order than is presently known for the anomalous dimensions 
of the twist- 2 operators which appear in the operator product expansion, which is fundamental 
to the mathematics underlying deep inelastic scattering. One reason that these are required 
is that they arise in the renormalization group equation for the evolution of the coefficient 
functions which are also required at a higher loop order. Indeed knowledge of the three loop 
operator dimensions and two loop coefficient functions will mean that the full evolution with 
energy scale can be performed at a new order. The anomalous dimensions have been studied 
at successive orders since the discovery that QCD was asymptotically free ||]. The original one 
loop calculations were carried out in |J whilst two loop corrections were constructed later both 
for the flavour non-singlet and singlet cases in @, |, g. These results were expressed not only 
as a function of the colour group Casimirs and number of quark flavours, Nf, but also of the 
moment of the operator denoted by n. To two loops the expressions as a function of n were 
reasonably straightforward to achieve. One requires this explicit dependence rather than, for 
example, a table of numerical values, in order to make contact with an alternative method of 
calculation. This is the splitting function approach of the DGLAP formalism, Q. The relation 
between the two being determined by the Mellin transform and its inverse where the conjugate 
variable to n is x which corresponds to the momentum fraction carried by the parton. Recently 
substantial progress has been made towards the determination of the three loop terms for the 
operator dimensions. In |Jsj, ^] exact expressions have been given for the even moments for the 
non-singlet, (n < 10), and singlet, (n < 8), operators. These calculations were performed by 
use of computer algebra and symbolic manipulation techniques to handle the huge number of 
Feynman diagrams which occur and to organise the resulting expressions for each graph. The 
corresponding coefficient functions are also known for the same values of n, j|, g|. Therefore 
it remains to construct the full n-dependent result. Clearly this is a highly non-trivial but 
important challenge. 

Some insight into the form of the operator dimensions has, however, already been provided 
through the large Nf technique, [H^, 11]. This is an alternative method to conventional perturba- 



tion theory where one can express all orders results in the coupling constant of a renormalization 
group function as a closed function of the spacetime dimension d. For QCD the method has 
extended the pioneering work for the O(N) a model developed in the series of articles ]jl2] , pj. 
Through ideas in the critical renormalization group one can directly relate the coefficients in 
the e-expansion of this function or critical exponent, where d = 4 — 2e, to the explicit pertur- 
bative coefficients at the same order in powers of 1/iVy. So, for example, using this approach 
an expression was derived for the /3-function of QCD, flij] , and the non-singlet twist-2 operator 
dimension at 0(1/Nf), |L0|]. The three loop coefficient of the latter, which was given explic- 
itly as a function of n, agrees exactly with the full calculation of j|] for the even moments 2 
< n < 10. Likewise the extension of that work to the singlet case, jLlJ], yielded similar exact 
agreement for the region of overlap, ||. Clearly these analytic results will play an important 
role for checking the full n-dependent three loop result. However, the singlet 1/Nf calculation 
of |ll[] only concentrated on one of the operators which occurs. In perturbation theory there 
is mixing between two operators which are respectively predominantly fermionic and gluonic in 



nature, y]. As the former was studied in [11] it remains to derive results for the latter type of 



2 



operator in order to complete the full leading order in 1/Nf computation. This is the main aim 
of this paper where we will discuss the technical issues required to obtain a result on a par with 
[11]. Although these are provided with the particular case of QCD in mind, the calculational 
algorithms to determine the values of the scalar integrals which arise are reasonably general 
and therefore applicable to the determination of the anomalous dimensions of similar singlet 
operators in other field theories in the 1/Nf method. Finally, we remark that although this 
paper is concerned with large Nf methods for the operator dimensions, some insight into the 



n-dependence of the coefficient functions in large Nf have already been provided in Jl^, 16]. 

The paper is organised as follows. We review the basic formalism and background to our 
problem in section 2. The details of the computation of the Feynman diagrams which occur in 
the large Nf critical point formalism are provided in sections 3 and 4. These deal respectively 
with the QED and non-QED sectors of the operator dimension defined essentially by those 
graphs which do not and do vanish by Furry's theorem. Section 5 is devoted to the derivation of 
explicit n-dependent results as well as several concluding remarks. The appendices respectively 
provide some basic integral results and an explicit three loop calculation to illustrate some of 
the internal cross checks which were necessary to validate the general recurrence relations we 
had to derive. 



2 Background. 

We begin by reviewing the necessary features of the perturbative structure of the anomalous 
dimensions of the twist- 2 singlet operators we are interested in. These are defined to be, 

pw-Mn = jn-l^ 7 W£)Wi . . . D^tp J - trace terms 

= \i n - 2 StiG a ^ u D^ 2 ...D^G a ^ -trace terms (2.1) 

where ip u is the quark field, AJJ is the gluon field, = 8^ + T a A*, G® v = d^A a u - d v A% + 
f a * 3C A b A^, T a are the colour group generators, f abc are the structure constants, 1 < i < N c , 1 < 
I < Nf, 1 < a < (N 2 — 1) and S denotes symmetrization in the Lorentz indices. In perturbation 
theory the canonical dimension of each of the operators is equivalent. Therefore they mix under 
renormalization which leads to a matrix of renormalization constants and associated matrix of 
anomalous dimensions. To fix notation we define the renormalization constants by 

Of n = Z l3 oY Te (2-2) 
where the anomalous dimensions, 7^(0), are defined by 

7«(a) = ( 7qq [ a \ M'l ) (2-3) 

with 7y(a) = (5{a){d / da) In Zij and (3{a) is the usual renormalization group function governing 
the running of the coupling constant a. The entries in 7«(a) depend on the colour group 
parameters, Nf and n and since it is the 1/Nf corrections which we will focus on we define the 
explicit form of each the entries as 

l qq {a) = aia + (ojjiJV/ + a 22 )a 2 + (a 31 Nf + a 32 N f + a 33 )a 3 + 0(a 4 ) 

lgq {a) = ha + (b 21 N f + b 22 )a 2 + {b zl Nf + b 32 N f + 6 33 )a 3 + 0(a 4 ) 

l qg {a) = c 1 Nfa + c 2 N f a 2 + (c 31 N f 2 + c 32 Nf + c 33 )a 3 + 0(a 4 ) 

lgg {a) = (d 11 A / + d 12 )a + (d 2 iA / + d 22 )a 2 + ( ( i 31 iV/ + ( i 3 2A / + ( i 33 )a 3 + 0(a 4 ) (2.4) 



3 



Clearly the Nf dependence does not enter in the same fashion in each term. The explicit 
determination was carried out in j|, || ^] but we record for completeness only those which we 
require 
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(2.5) 



where the colour group Casimirs are defined by T a T a = C 2 (R), tr(T a T b ) = T(R)5 ab and f acd f bcd 
= C 2 (G)5 ab . The function Si(n) is defined by Si(n) = J2r=i l/ r '- These expressions ( fE5| ) are 
determined in perturbation theory by inserting both operators, fl2.l|), in quark and gluon 2-point 
functions and deducing the pole with respect to the dimensional regulator before using MS to 
define Z%j. The procedure to find leading order 1/iVf information on 7ij(a) is similar and has 
been discussed in [10, 11]. Briefly one uses critical point renormalization group techniques to 
analyse the scaling behaviour of the operators we are interested in at a fixed point in the (5- 
function of QCD. This is defined to be a non-trivial zero of the d-dimensional /^-function and 
the value of the coupling is denoted by a c . As has been observed in [14] the propagators take a 
particularly simple scaling form at a c . Further, the anomalous dimension of each operator can be 
deduced in powers of 1/Nf by noting that there is a relationship between the critical exponent of 
the operator and the analogous anomalous dimension renormalization group function evaluated 
at a c . In the (i-dimensional case such a non-trivial zero exists and, moreover, is computable in 



a power series in 1/Nf. For instance, from the 3-loop result, |17|, 
(d-4) 



], in the notation of [|19| , 



(3(a) 

+ 

one deduces that as Nf 
3e 



-a + 



~T(R)N f - l ^C 2 {G) 



20 34 
4C 2 (R)T(R)N f + -C 2 (G)T(R)N f - -C 2 2 (G) 



a 3 + 0(a 4 



(2.6) 
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(2.7) 



where we have included the contribution from the recently computed exact 4-loop MS /3-function, 
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The scaling behaviour of the quark and gluon propagators in momentum space at a c are 



then simply, [ 10 1 



- (1-6) 



2i 



where 6 is the covariant gauge parameter, d = 2/j, and the exponents a and (5 are defined as 

a = n - 1 + ±7? , = 1 - v - X (2-9) 

The quantity 7/ is the quark anomalous dimension which corresponds to the quark wave function 
renormalization at the critical coupling and x is the- anomalous dimension of the quark gluon 
vertex. The canonical dimensions are deduced by demanding that the d-dimensional action is 
dimensionless. This is contrary to perturbation theory where the canonical part is defined by 
ensuring the 4-dimensional action is dimensionless. The quantities A and B are momentum 
independent amplitudes. Using (|2.8| ) various leading order results have been established such 
as the QCD /3-function at 0(1/Nf), []l4|. This encodes the coefficients of all powers of e in the 
0(1/ Nj) term of (|2.7| ). As the leading order values of rj, x an d the amplitude combination z = 
A 2 B are required here, we note their explicit values are, pfl], 

[(2/j-l)(/i-2)+/z6]C 2 (i?)7ff 
Vl = —r, l 2 - iu J 



(2/i-l)( M -2)T(i2) 
-!)(//- 2) +/xfr]C 

(2/i-l)(M-2)T(i?) 2(2 M -l)(/x-2)T(i?) 

Of + M 
2(2 M -l)(//-2)T(i2) 

where, for example, 77 = Ya^i Vii 6 )/^* an d 



[(2/x -!)(//- 2) +//&]C 2 (fl>ff [(2»-l)+b(»-l)]C 2 (G)rg 

Xl = W -DO" ^ ^ ■ ^ — (2 ' U) 

* = , r( ^ +1)?? L^ (2-12) 



o _ (2/x-l)(/i-2)r(2/i) 

771 " 4r^)i> + i)r(2-^ 



In [11 1 it was these propagators which were used to analyse the fermionic operator O q and 
produce an expression which agreed with known results to three loops. This remark needs to be 
qualified, however, since the canonical dimensions of the fields at a c differ from the perturbative 
ones. Therefore the canonical dimensions of ( |2.1[ ) are different and so they cannot mix under 
renormalization. In other words the mixing matrix at criticality is triangular at least at leading 
order in 1/Nf and both operators cease being twist-2 at a c . Instead O g is twist- (2 + 0{e)). 
Although this would appear to mean that the critical point large Nt method fails to have any 
contact with perturbation theory which we are seeking, there is a relation between 7ij(a c ) and 
the scaling dimensions of (2.1) at a c . This was established in []ll[] where it was noted that the 



eigen- anomalous dimensions of (2.3) evaluated at a c are in fact the dimensions of the operators 



(2.1) whose fields have the canonical dimensions given by (|2.9| ). In other words in perturbation 



theory the information contained in the critical exponents relating to the dimensions of (2.1) at 
a c are given by two combinations of the perturbative coefficients listed in (|2.5| ). To make these 
comments more transparent we return to Jij(a) and deduce the eigen-anomalous dimensions at 
a c in powers of 1/Nf. First, the eigenvalues for general a are given by 

1 1 r if 

A±(a) = -{iqq + lgg) ± - [(jgq ~ Iggf + ^Hqglgq] ' (2.14) 
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However, expanding in powers of a and retaining the same orders in 1/Nt with the definitions 
(O), we find 
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AfV + 0(Nfa 



(2.15) 



Noting that a c is 0(1/Nf), it is easy to observe that A_(a c ) = 0(1/Nf) whereas A + (a c ) = O(e) 
which bears out our remarks on the difference in twist at a c . As the dimension of the pre- 
dominantly fermionic operator has been discussed elsewhere, we record that the predominantly 
gluonic eigen-operator dimension is 



A+(a c ) 



3dn£ 



1 T3 / , &1C1 11 , , , 

+ 4 f rf 2i + ^ - (3C 2 (fl) + 5C 2 (G)) du ) e 2 
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+ o 



(2.16) 



We note that the location of a c is present in this expression because of a non-zero du term at 
one loop in (p 



Two technical issues remain to be discussed which are fundamental to the calculation of 
A + (a c ). The first relates to a property of the fixed point. At a c it was demonstrated by Hasen- 



fratz and Hasenfratz, |21], that at 0(1/Nf) QCD is equivalent to a simpler lower dimensional 
model which is the non-abelian Thirring model, (NATM) . This means that they lie in the same 
universality class and therefore critical exponents are the same in each model. In other words it 
does not matter which theory we study at a c , the same d-dimensional expressions are obtained 
at leading order. Similar and more widely studied d-dimensional equivalences exist like, for 
example, that between the O(N) a model and O(N) 
completeness the respective lagrangians are 



theory. (See, for example, p2|.) For 



L QCD 



L 



NATM 



4e 2 
2A 2 



(2.17) 
(2.18) 



where e and A are the coupling constants which are dimensionless in 4 and 2 dimensions re- 
spectively. However, since the NATM has fewer interactions due to the absence of cubic and 
quartic gluon self-interactions it is more practical to use it. The effects of such interactions are 
recovered in diagrams with quark loops with respectively three and four gluon legs, pl"[ . Indeed 
this equivalence has already proved useful in establishing several results in large A^- QCD, 

The second point concerns the regularization of the theory at the fixed point. Unlike per- 
turbation theory we use a fixed non-regularizing dimension d in the critical point approach and 
therefore the only sensible type of regulator which remains to us is an analytic one. This is 
introduced by shifting the gluon dimension (3 to (3 — A where A is a dimensionless infinitesimal 
parameter, [13|. It plays a similar role to the e-regularizing parameter of dimensional regu- 



larization. Therefore when we insert O g in a gluon 2-point function and replace the lines by 
the critical propagators with a non-zero A the value of the diagram is no longer infinite but 
proportional to 1/A. As has been discussed in [| 



11 



the residue of this simple pole is an n 
and d-dependent function and it is this which contributes to the full value of A+(a c ). 
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3 QED sector. 



We now turn to the discussion of the technical issues involved in computing the critical exponent 
A + (a c ) beginning with some general considerations. As in a perturbative calculation we insert 
the composite operator O g into a gluonic 2-point Green's function and compute its divergent part 
with respect to the critical regularization. This is illustrated in fig. 1 and the total value of such 
contributions correspond to the renormalization of the bare operator whose critical exponent is 
rjOg- To obtain the full value one must also include the wave function renormalization of the 
fields present in the composite operator as in perturbation theory. Since these are gluons then 
the dimension is 

A+(a c ) = - V ~ X + Vo„ (3.1) 



It is worth noting that because the operators ( |2.l[) are physical the combination (3A) is inde- 
pendent of the gauge choice even though each of the individual quantities of ( |3.1[ ) depend on 
b. Therefore we have chosen to perform the calculation in the Feynman gauge. From [^, ||] and 
(|2.1| ) the Feynman rule for the (bare) insertion which is quadratic in the gluon field, is, for an 
external momentum p, 



l + (-l) n ] 7]^(A P ) 2 - {A^ Pu + A uPll )Ap + p 2 A IM A v (Ap) 



(3.2) 



where /U and v are the free indices which remain after contracting with the null vector A^ 
of the light cone projection, [11]. (Although we have used the same symbol for the analytic 
regularization there ought not to be any confusion in what follows as one is clearly a scalar 
object whilst the other is a Lorentz vector with the property A 2 = 0.) The Feynman rule for 
the insertion with three or more gluon legs have been given elsewhere, Q. At the order in 1/Nf 
which we are considering there are only a few contributing graphs whch are illustrated in fig. 2. 
There are no graphs with self energy insertions on internal lines since we are computing with 
propagators which have non-zero anomalous dimensions. Each of these graphs in fig. 2 will in 
principle have the same general Lorentz structure in its pole part which is 



A Vflu (Ap) 2 + B(A fl p u + A uPf ,)Ap + C^(Ap) 2 + Dp 2 A^A 



(Ap) n - 



(3.3) 



where A, B, C and D will depend on /i, n and the regulator though it is only the residue of the 
simple pole in the latter which we are interested in. In principle we ought to compute each of 
these expressions for each of the integrals, which would require taking various projections on the 
explicit form of each graph. However we can minimize the amount of computation that needs 
to be performed by appealing to the renormalizability of QCD. As the theory is renormalizable 
then the sum of the contribution of the form ( |3,3j ) from each of the graphs of fig. 2 must be 
proportional to (|3.2| ) . The resulting coefficient of proportionality is then simply related to rjo g ■ 
Therefore it is much simpler to compute only one of the four quantities in ( |3.3| ). (However, we 
note that we have checked this assumption by explicit computation for the QED sector of fig. 2 
which involves only the first two 2-loop graphs.) We have chosen to deduce A and it is a simple 



exercise to invert (|3.3|) to produce 



A 



2(M-1) 



^{ApY - {A^pv + A uPfl )Ap + p 1 A il A v {Ap) n - 2 Y 



fJ,U 



(3.4) 



This projection on each of the four graphs can be simplified further by noting that renormal- 
izability also requires that C is finite with respect to the regularization. In the same inversion 
which produces (|3.4j) . 



C 



2(/i - 1) 



(3.5) 
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Therefore the projection of each of the integrals we take to produce the divergent part with 
respect to the regularization is 

[ V^(Ap) 2 - (A„p v + A uPfl )A P ] (Ap) n ~ 2 rV (3.6) 



2(/i-l 

This step substantially reduces the number of constituent integrals which would have to be 
considered in, for example, the final graph of fig. 2. A second simplification of the organization 
of the calculation is to examine the contributions from the sectors defined by the colour group 
Casimirs C2{R) and 02(G) separately. As the latter arises due to non-zero colour group structure 
constants we will refer to the former as the QED sector and concentrate on determining its 
contribution first. One reason for this is that there are various technical issues in the full 
calculation which we need to discuss and these are best illustrated in this case. In the 02(G) 
sector other techniques will be introduced which will assume these results. 

For the remainder of the section we will discuss the first two 2-loop graphs of fig. 2 as only 
they involve C2(R)- The first of these graphs is elementary to determine, primarily because it is 
equivalent to chain integrals in the language of However the presence now of an operator 
with a non-zero moment means that more general chain results are required. These have been 



given in appendix A as they have been widely used by other authors in this area, [23]. Therefore 
we merely record the total contribution from this graph to rjo g is 

[fi 2 n 2 + /i 2 n - 2fi 2 - 9/in 2 - 9/rn + 2/j, - n 4 - 2n 3 
+ 9n 2 + 10n]n(jJi - l)T(ji)T(n + 2- n)C 2 {R)rR 

/[n(n + l)(p -2)(n-n- l)r(3 - n)Y(ji + n + l)T(R)} (3.7) 



The computation of the second integral of fig. 2 is more tedious to determine. From Q3.6| ) 
we treat the contributions from the rf v and (A^p u + A v p^) projections separately. To proceed 
one contracts with either projection tensor and then takes the trace over the quark loop which 
leaves a set of two loop scalar integrals. The majority of these integrals can be calculated in the 
same way as those of the first graph of fig. 1 since they are elementary chains. However a set of 
two loop graphs remain which have the form illustrated in fig. 3 which clearly require another 



technique. We use a coordinate space notation for integrals, similar to [12], but the explicit form 
of that graph is given by 

r t n \ - ( (AyY(AznA(y-z)Y 
mn[ ] ~ J yz y 2 (x - y)W(x - zf((y - *)=»)« ^ 

where p, q and n are integers and a is an arbitrary exponent. Here and elsewhere we will use 
the notation that a line in a graph which has a bracketed integer, (p), beside it corresponds to a 
factor in the numerator of p contractions of its vector with A„. Unbracketed symbols beside a 
line correspond to its propagator exponent. From the nature of the original integral the central 
(y — z) line corresponds to the location of the operator insertion with moment n. It turns out 
that only several values for p and q are needed for n arbitrary. These are (p,q) = (0,0), (1,0), 
(1, 1) and (2,0). However, some of these are related by symmetry or by rewriting Ay as Ay = 
A(y — z) + Az. Therefore for n even 

Gion(a) = ^G 00 n(a) (3.9) 

and for arbitrary n 

Gun = G 2 on(a) - ^G 00 ( n+2 )(a) (3.10) 

It therefore remains to compute Goon(a) and G20n(ct) which is achieved by establishing a recur- 
rence relation for the more general integrals where the contractions of all numerator vectors 
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has not been performed. Then we will solve this for the cases we are interested in. It is best 
to illustrate this for the simpler Goo n (a) case first. Then using Lorentz symmetry that integral 
will decompose into the form 

A(a,n)x fll ...x^ + B(a,n)[r) IJrllX2 x IJl3 ...x IMl + {\n{n - 1) - 1) terms] + 0(x n ~ 4 ) 

(3-11) 

In this and similar decompositions we will omit the overall power of the factor of (x 2 ) as it plays 
no significant role in the derivation and can readily be restored by simple dimensional analysis. 
Here A(a, n) and B(a,n) are invariant amplitues and our aim is to determine a recurrence 
relation involving A(a, n) only. This is because after contracting (|3.11| ) with A'' 1 . . . A^™ only 
A(a,n) remains which is the value of Goo n (a). The amplitude B(a,n) is necessary to achieve 
this as can be seen in the explicit derivation, though the terms in powers of x will not be 
necessary and these have been indicated by the order symbol. In the following construction of a 
relation for A(a, n) we consider only the leading contribution in powers of x when we take two 
contractions. These are r/' 1 " -1 ' 1 " and x^ n and give respectively 

A(a-l,n-2) = A(a,n) + 2(fi + n - 2)B(a, n) 

2Goog-i) = A(a,n) + (n-l)B(a,n) (3.12) 

where G~^ (a) corresponds to the integral of fig. 3 but with the exponent of the top left line 
replaced by zero. In other words it is an elementary chain integral of the form which has arisen 
already. Prom these two expressions we deduce, for arbitrary a, 

G o(n-2)(a - 1) = — ( n - 1) — G oo(n-i)(«) / _ -jn — Goon(a) (3.13) 

from the decomposition the various integrals we require are Croon(2/^ — p + A) for p = 2, 3 and 
4. Indeed once we have obtained one of these values the others will follow from elementary 
application of ( p.l3| ). It turns out that Goon(2/U — 3 + A) is A-finite. The justification for this 
will prove to be fairly general and will be used later in solving recurrence relations for other 
integrals. Essentially it follows from the fact that the dimension of the integral is greater than 
unity. We first recall the usual calculational technique used in [^] for computing integrals 
without a moment factor in their numerators but which are divergent and similar in form to fig. 
3. The first step is to transform it to the momentum representation defined in In this one 
has to include factors arising from the Fourier transform. In particular one factor is related to 
the dimension of the integral. If this has the form — [p + A) for a positive integer p, then the 
corresponding factor is divergent with respect to A. It therefore remains only to compute the 
momentum representation integral for A = 0. This is possible as the final form is related, for 
example by integration by parts, to known integrals like those of Chetyrkin and Tkachov, [p4|| , 
denoted by ChT(a, (3), [Oj. It follows therefore for this case that if the original integral has 
a dimension greater than unity then it is finite. For the integrals with a non-zero moment the 
argument is similar. Indeed we checked this reasoning by explicitly computing Goon(2/i — 3 + A) 
for various values of low n. For future reference we make some comment on the extension of 
this argument to other integrals not necessarily of the exponent structure of fig. 1. It may not 
always be the case that the corresponding n = integral is A-finite even though its dimension 
satisfies the above criterion. However, its moment partner may be. This is because the presence 
of a numerator factor like A(y — z) alters the nature of the singularity which we are interested in. 
These remarks can be simply illustrated by considering the elementary function a n (a) which is 
common in all our calculations. For example, the function oo(/x — A) is divergent in the limit A 
— ► 0, whereas for a non-zero n a n (/x — A) is finite in the same limit. The presence of the non-zero 
moment appears to be regularizing the infinity of the n = function and it is this which occurs 
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in the two and three loop integrals. Indeed throughout our application of recurrence relations 
we always have used this dimensional argument for integrals to determine the terminating 

case. For all the cases we are interested in we checked the finiteness explicitly for at least one 
value of non-zero n. The argument appears robust. 

Therefore given this property we can deduce Goon(2/i — 4 + A), for example, from ( |3.13|) 
which now is reduced to a computation of chain integrals. Again we have checked this result by 
carrying out the explicit calculation for low n and found agreement. To complete the calculation 
of the second graph a recurrence relation for G20n(«) is deduced in the same way as for (gUp . 
The major difference here is that the decomposition analogous to (3.11) now has four terms 



to 0(x^ n ~ 2 ^). However there are various contractions which can be made from which one can 
deduce, for example, G ! 20n( Q! )- Finally the two projections of ( |3.6| ) can be performed and the 
contributing value of the second integral to rjo g is 

- [8fi 4 n 2 + 8/i 4 n + 8// + // 3 n 4 + 2//n 3 - ll//n 2 

- 12// 3 n - 24^ 3 - 3fi 2 n 4 - 6^V - 19^V - 16^ 2 n 
+ 24/i 2 — /un 6 — 3^n 5 — 5/in 4 — 5/un 3 + 34/m 2 + 36/m 
-8fi + 2n 6 + 6n 5 + 10n 4 + 10n 3 - 12n 2 - 16n] 
x M rOu - l)r(n + 1 - n)[C 2 (R) - C 2 (G)/2]7?° 

l\n(n 2 - l)(n + 2)r(3 - n)r(p + n + l)T(R)} (3.14) 
We close this section with the remark that we have made use of the computer algebra and 



symbolic manipulation packages Reduce, |25|] , and Form, [|26| , in this and subsequent graph 



evaluations as there is a substantial amount of tedious algebra to handle. 



4 C 2 {G) sector. 

We now detail the calculation of the remaining two graphs of fig. 2 each of which involves only the 
colour group factor C*2(G). The strategy is the same as in the C*2(i?) sector. In other words we 
take the trace over the fermion loop and evaluate the pole part of the resulting bosonic integrals. 
The major effort in the calculation resides in the determination of the integrals which are not 
reducible to chains. First, we consider the two loop graph. The three gluon insertion Feynman 
rules can be found in, for example, 0|. It has a part which involves a sum over the operator 
moment and it is therefore not clear whether it will be possible to obtain a closed analytic 
expression for the contribution from this graph. Moreover, the distribution of the exponents 
around the two loop topology is different from that of fig. 3 and so a new recurrence relation will 
be needed. The integral is defined in fig. 4 where a and j3 are arbitrary exponents and n and p 
are integers. In addition to the simple chain diagrams given by the graph of fig. 4 where each 
of the lines with a unit exponent is successively replaced by zero, we have to deal with those 
graphs where the exponents of each of the lines are also replaced by an exponent of (— 1). These 
integrals are easily related to a sum of chain diagrams by integration rules like that illustrated 
in fig. 5 which is simple to establish. As these resulting chain integrals are elementary to deal 
with we focus on the construction of a relation for H pn (a,(3) and consider Ho n (a,f3) first. A 
clue on how to proceed comes from the case Hoo(a, j3) which has been determined in closed form 



in [24|. It is related to C/iT(/x — a,/i — (5) of [|12J after applying the momentum representation 
transformation of p2j| . The presence of a non-zero moment n however, disrupts the construction 
of a closed form when the same algorithm is applied. Instead we apply a series of transformations 
to HQ n (a,f3). Noting that under a conformal transformation — » j z 2 and therefore Az — > 
Az/z 2 , we apply a left conformal transformation on fig. 4 with the origin at the left external 
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point. The lines joining the upper vertex now have no moment dependence and so we can apply 

a 



transformation of [12]. Finally a momentum representation transformation leaves us with 
an integral proportional to the graph given in fig. 6. Again as the lower vertex has no lines 
joining it which have a moment dependence we can integrate by parts on it using, for example, 



the rule of [23] and relate it to four integrals. Three of these are simple chains whilst the fourth 
is related to Ho n (a,(3 — 1) when one undoes the original transformations. The result therefore 
is, for arbitrary a and (3, 



H 0n (a,P) 



(2/i - a - P + n - l)(a + /3 - (i) 



H 0n (a,P-l) 



09- l)( M -/9-l) 
a(l)o(/3 + 2 - (i)a n {3(i -a-(3-2 + n) 
{(3 - l)a(3 - /i)a„(2/u - a - 2 + n) 
09-1 



(/X-/3-1) 



i/(2, 13,2(1-/3-2) + v non (a,2,2(i-a-2 + n) 



(4.1) 



For ( |4.1| ) to be practical we once again need to have a terminating value which will be an 
integral which is finite with respect to the regularization. Unlike in section 3 we require several 
distinct series of expressions such as Ho n ((i — 1 + A, (i — 2 + A), H$ n (2(i — 3 + A, A) and 
Ho n (2(i — 3 + A, (i — 1 + A). For the former the natural endpoint is Ho n ((i — 1 + A, (i — 1 + A) 
= O(l). This has been verified for various values of small n and simple substitution reveals 



flbn(M- 1 + A,M-2 + A) 



(2/x-3)(//-l) t/(l,l,2//-2) 
2(ra + l)(/x-l + n) r(/x)A 



(4.2) 



which we have also checked for low n. Of course we have only recorded the value of the 0(1/ A) 
part of the integral since the contribution to rjo g only arises from the residue of the simple pole. 
One cannot always directly use (4.1) for each integral due to the fact that if, for instance, (3 = 
(i — 1 + A then a denominator factor of (4.1) is A. This would require the associated two loop 
integral to be computed to the next order in A which would be tedious. Instead one can employ 
the identity 



H 0n ((L-2 + A,/x- 1 + A) = 



n 



(-l) r ff0r(M- 1 + A, (JL-2 + A) 



(4.3) 



r=0 



and insert the previous result. Elementary algebraic identities allow the finite sum to be deter- 
mined, leaving 



H 0n ((i-2 + A,(L-l + A) 



(2/^-3)^-1)1/(1,1,2^-2) 
2((x - 2)r(/i)A 



1 



r(M-i)r(n + i) 



(n + l) 



T((jL + n) 



(4-4) 

For the other structures we note that by our earlier dimensional argument Ho n (2(i— 3+ A, 1 + A) 
and Ho n (2(i — 2 + A, A) are A-finite and lead to the following results 



H 0n (2(i - 3 + A, A) 
Ho n {2(i - 4 + A, 1 + A) 
H 0n (2/j-2 + A,-l + A) 



i/ 00n (l,2,2/i - 3) 
2(n + l)T((i + n)A 

(2// - 3)i/(l, 1,2/x- 2) 
2((i-2){n + l)T((i + n)A 



r(4 - n + n) 

r(3- M ) 



T(n + 2) 



(/x- l)(2/i-3)i/pon(l,2,2/z-2) 
2(n + l)r(A* + n)A 



(4.5) 



Again we have verified these (and others) all explicitly for at least one non-zero value of n. It 
turns out that the remaining integrals which we need of the form Ho n (2(i — p + A, (i — q + A) 
for positive integers p and q are A-finite. 
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To complete the construction of results for the topology of fig. 4 we note that applying the 
generalized moment integration by parts rule of |23| to -£f( p -i)( n +i)(a, /3) gives the reduction 
formula, for arbitrary a and 0, 



H pn (a,(3) 



(2/x-a-/3 + n-l) 



(n + 1) 



■«(p-l)(n+l) (a, (3) 



a 



(n+1) 

P 

(n+1) 



H (p-l)(n+l)( a + 1) P) ~ H (p-l)(n+l)( a + ^) 
H^ 1){n+1) (a,P + 1) - ff^ lKn+1) (a,/3 + 1) 



(4.6) 



where the superscripts on an //-function here refer to the exponents respectively of the top left, 
central and top right lines relative to the diagram of fig. 4. The result ( |4.6| ) allows us to reduce 
any integral H pn (a,(3) with p > down to ones of the form Ho n (a,(3) whose values are known 
from ( [4,ip . Therefore one can now assemble all the pieces to both projections of the original 
two loop diagram. 

For that part of the operator insertion which does not involve a finite sum over moments 
one obtains the contribution 

(4//n 2 - 4// - 10^ 2 n 2 + 2/u 2 n + 8u 2 + 8nn 2 - Ann -3a- 2n 2 + 2n)T(a + n) 

- (2u 2 n - 2u 2 -2>un + Au + n- l)T(u + l)T(n + 1)1 rftC 2 (G) 
/[2(2u - - 1)0* - 2)(n - l)T(u + n)n] (4.7) 

However, the total of all pieces to that part which does involve the sum yields, ignoring colour 
factors and powers of z for the moment, 



4(/x - l)i/(l,l,2/x - 2 



n-2 

■E 



(2u - l)r( M ) ^ r(r + l)r(/i + r + 1) 
x "[4r(/i - l) 2 + l]T(n + r + 1) - [r(2a - l)(/x - 1) + /x)]r(/x)r(r + 2)] 

From the definition for 5j(n) and the result 



A £(r + a) 

^ r( 



l 



r=l 



1) 



r(a + i) 
r(6) 



T(n + a + l) 



is 



(2/i-3)(2/x-l) Si(n) 



\r + b) (a-b 

then we deduce the closed form for 

4(/x - 1,2/i- 2) 

(2/i - i)r(At) 

a / x (2n - 1) 
n(n — lj 

Thus the total contribution from this two loop graph is 
2n(n-l)Sx(n) (2un-n + l)T(u + l)T(n) 



T(n + b) 



(4.8) 



(4.9) 



(n + 1) 



n 



2+[(2 M -l)(n-l) + l] 



r(n-i)r(^) 

T(a + n-l) 



(2u-l)(u 



+ 



(2n 2 n — 3/un + \i + 2n) 
2(2/i - l)(a - l)(u - 2)T(n + n) ~ 2(2u - l)(u - l)(u - 2)n 



(4.10) 



rRC 2 {G) 
T(R) 
(4.11) 



To complete our computation we now turn to the final graph of fig. 2, which is the hardest of 
the four to evaluate. As we are in effect performing two separate projections on the integral itself, 
we note that multiplying the graph by A^p u , for instance, allows us to immediately complete 
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the integration for the quark loop that the p v contracts into. The resulting two loop integral is 
effectively completed since its constituent scalar integrals have already been evaluated for the 
third graph of fig. 2. A similar property arises for the rj^ u contraction of the second term of 
the operator insertion (|3,2| ). As for the A^p u contraction there is also no new feature for this 
case. Instead we consider the remaining two terms of Q3,2j ) in fig. 2 and their rj^ u projection. 
In taking the relevant traces one obtains scalar integrals which are now either simple chains or 
two loop integrals which have already been considered. In addition we obtain a new three loop 
topology whose most general form is given in fig. 7 and has a similar integral definition to (|3.8| ). 
Here p, q, n and 7 are integers with a and (3 arbitrary exponents. We emphasise that here the 
line with exponent ( — 7) corresponds to a numerator factor of [(y — z) 2 ] 7 . Similar to the other 
02(G) two loop chain type graphs there are a variety of three loop graphs of the form of fig. 7 
but where one of the lines with unit exponent is replaced by ( — 1) for 7 = 1, 2 or 3. These can 
readily be reduced to at most two loop graphs already computed using integration rules similar 
to that of fig. 5. 

Therefore since the structure of fig. 7 is new to the calculation we give some details of its 
treatment, though the algorithm is no different from those already established for the two loop 
cases. We consider the construction of a relation to deduce Coo 7 (ci!, P, n) first and note that its 
decomposition into Lorentz invariants where the vectors have been omitted is the same as in 



( 3.11 ). However, due to the form of the present integral the contractions •q^-i^n anc [ x v-n yi e ld 
respectively 

A(a-l,p,n-2) = A(a,/3,n) + 2(ij, + n-2)B(a,p,n) (4.12) 
A(a,P,n) + (n-l)B{a,P,n) = \[A{a, p, n-1) + A(a-1, p, n-1) - A(a, (3-1, n-1)] (4.13) 
Eliminating B(a,(3,n) gives 

A(a-l,P,n-2) = - {2fl + n ~ 3) A(a,f3,n) 

(n-1) 

+ V n ~ 2 ' [A(a,f3,n-l) + A(a - 1, p, n - 1) - A(a,(3 - l,n - 1)] 
(n-1) 

(4.14) 

This is not quite in a form which is practical since we require a relation between A(a, f3, n) and 
say A(a — 1,(3, n) where the former set can be set to zero upon reaching a terminating value. 
The presence of the final two terms in the square brackets leaves us an equation which will relate 
several integrals without allowing us to evaluate them individually. To eliminate these problem 
terms we consider the general integral of fig. 7 and apply a conformal left transformation to 



it in the language of [12]. Then integrating by parts on the central 4- vertex with the line with 
exponent [3 as reference and applying the inverse conformal left transformation to the resulting 
five integrals, we obtain the general result, for arbitrary a and (3, 

C pqi {a,(3,n) = — — j— ^— ^^-[C pqy (a, (3 - l,n) - C pqy (a - 1, f3,ri)] 

' -[0«' P-l,n)- C™(a - 1,(3, n)} 



Ti 

+ ( a _ pj C pq7 (a-l,p,n-l) (4.15) 

where a similar notation to ( |4.6| ) is used but here the superscripts refer only to the exponents 
of the lines joining the top vertex. Although we are dealing with the leading term A(a, P, n) 
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in ( 4.14| ) we can use (|4.15| ) to eliminate the unnecessary terms since after the contraction with 
A^ 1 . . . A Mn only A(a,f3,n) will survive and be equal in value to Cooy(a, P,n). Consequently 
after some elementary algebra which includes shifting exponents, we deduce, for arbitrary a and 

P, 



Coo 7 (a,/3,n) 



{2n + n - l)(2n - a - (3 + n - 2) 

(n + l)(/i-a-/?-2) 
(fi + n)(2fj,-2a + n-3) 



C 007 (a + l,/3,n + 2) 



(n + l)(p-a-/3-2) 

2(/i + ") r^Gi 



C o 7 (a + l,P,n + 1) 



(n +!)(//- a 2) 



[C u ^(a + 1, /? - 1, n + 1) - Cgg; (a, /3, n + 1)] 



(4.16) 



This differs in a minor respect from the two loop relation of ( |4.1| ) in that the first two terms on 
the right side are not at the same moment level. However, the dimension argument we use to 
find a terminating value is valid for all n > 1 and these terms will be A-finite. 

Relations are also required for C\o 1 {a, (3, n), Cn 7 (a,/3, re) and e2o 7 (a, /3, n) and these are 
determined by the same method as for Coo 7 («,/?, re) though there are of course more invariant 
amplitudes which can appear. However, one constructs the result for Cio 7 (a,/3, re) first since 
the decomposition of the other two relies on this lower level expression and Coo 7 (ck, P, re). For 
completeness we quote the recurrence relations for each level in appendix A where again the 
integrals where a line is zero are easily computed. As an aid to an interested reader we record 
the explicit results of the recurrence relations for a set of representative integrals at each level, 
noting that we have used the obvious values such as Cqq${2u — 3 + A, a — 1 + A, re) = 0(1) for 
all re > 1 to terminate the relations. We found 



O003 (2/x — 2 + A, /i — 1 + A, n) 



2(2 / Li-3)( / u 2 -l)rV + re + l) 

-T(n + 2)T(n + 3)]a(2/x- 2) 

+ 2(2/i - 3) (a* - l)a n (2/^ - 2)r(p + 2) 

x a(2u - 2)a^{l)u/[A{2u - \) 2 {u - 2)] 
/[(re + 2)(re + l)T(p + re)I»A] 



(4.17) 



Cioo(2/i - 4 + A, ix - 2 + A,n) 



(u 2 + 2,ure 2 + Sun + 4/i - 3re 2 - 12re - g)(ji - n - 2) 

x (/i-n-3)a„(2/i-2)r(u) 

+ (ji- 1) 2 (m - 2)(n + 2)a(2^ - 2)r(yu + re + 1) 

x (2/i - 3)a(2/x - 2)r 3 ( M )/[2( M - l) 4 (/x - 2) 2 ] 

/[(re + 3)(n + 2)(re + l)r(/i + n + 1)A] (4.18) 



C n2 (2fj, - 3 + A, fj, - 1 + A, n) 



\{2n 4 n l + 10/Tn + KkT 



o 2 2 

3// re 



19^ 2 re - 2bu 2 



an 



un + 7u + re 2 + 7re + 6)r(/Lt + n + 3) 



(/i 2 n + 2/^n 2 + 4/in + a 



re 



3re - 2)(/i - 1) 



x r(//)r(n + 5)]a(2/i - 2) 
+ (/i 4 + 2,u 3 re 2 + 8^i 3 n + 4/i 3 + 4//V + 36//V 
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+ 110/iV + 132/i 2 n + 53/i 2 



4/m 4 - 38/m 3 



- 124/m 2 - 162 /j,n - 70/i + n 4 + 10n 3 + 35n 2 
+ 50n + 24)(2/i - 3)(/x - n - 2)a n (2/t - 2)r(/x) 

x a(2 M - 2)r 3 ( / u)^/[4(2 M - l) 2 (/i - l) 4 (/i - 2)] 

/[(n + 4)(n + 3)(n + 2)(n + l)r(/i + n + 2)A] (4.19) 



and 



C201 (2/i — 4 + A, /i — 1 + A, n) 



[(/i 3 + 4/x 2 n + 5/i 2 + 5/ure 2 + 13/m + 8/z - 2n 2 

- 6n - 4)(n + 4)(n + 3)r( / u)r(n + 2) 

- 2(jm + 3/Lt - l)r(/i + n + 3)/i](n + 2)a(2/i - 2) 
- 4(/i 3 + 2/i 2 n 2 + 10/i 2 n + 9/i 2 + 2/m 4 + 20/tn 3 

+ 69/m 2 + 95/m + 44/x - n 4 - 10n 3 - 35n 2 - 50n 



- 24)(/x - n - 2)(/i - n - 3)a n (2/i - 2)r(/i) 

x a(2/i - 2)r 3 (/i - l)/[4(2/x - 1)(m - 2)(n + 4)] 
/[(n + 3)(n + 2)(n + l)T(fi + n + 2)A] 



(4.20) 



With the results from these recurrence relations and the earlier two loop results we summarise 
this intermediate step in determining the value of the three loop graph by giving the contributions 



from the r]^ u projections of the first and third terms of the operator insertion, (3^2). These are 
respectively 



[16/i 6 



16/x 5 n 



104/x 5 



4//n 2 + 60//n + 264// + 16// 3 n 2 



80//n 



+ 51/i 2 n + 232/i 2 + 11/m 2 



20/m - 88/z 



x W °C 2 (G)/[8(2/, -l) 2 ^ 



336/i 

- 2n 2 + 4n + 16]r(n + 2 - /x)r(/i - 1) 
2) 2 (n + 2)(n + l)r(2 - /x))r(/i + n)T(J?)] 



21/rV 



(4 M - l)r( M - l)r(n + l)ffl?C 2 (G)/[4(2^ - 1)( M - 2)r(/x + n)T(i?)] 
[16/i 6 + 16/i 5 n - 88/i 5 



4/x 4 n 2 



84//n + 192// + 12//n 2 + lU^n 



216// 



149//n + 136// + 8/m 2 + 63/m 



48// - 4n 2 - lOn + 8]/try?C 2 (G) 



/[8(2/i - l)(/x + n - - l) 3 (/i - 2)(n + 2)(n + 1)T(J2)] 



13/iV 



(4.21) 



and 



160//n 3 



+ [16//n 2 + 48// n + 32// + 32//n 3 - 56//n 2 - 360// n - 256// + 16//n 4 

■ 5 88//n 4 + 312//n 3 + 614//n 2 - 1518//n - 1530// 

.3 ioo..2„4 i ino..2„3 



- 72// n 2 + 1048// n + 852// 

+ 184//n 4 - 288//n 3 - 1145//n 2 + 1141//n + 1600// - 182//n 4 + 108//// 
+ 973/i 2 n 2 - 401/i 2 n - 970/i 2 + 83/m 4 + 4/m 3 - 385/m 2 + 34/m + 312/i - 14n 4 

- 8n 3 + 54n 2 + 8n - 40]r(n + 1 - n)r(p - l)/ir/°C 2 (G) 
/[8(2/x -l)(n- 1) 2 (m - 2) 2 (n + 2)(n 2 - l)r(2 - n)T(p + n)nT(R)} 
(2/i 2 + 3/m - 3/i - n + l)r(/i - l)r(n - l)nrftC 2 (G)/[2(2/x - l)(/i - 2)r(/i + n)T(R)\ 
[16// 7/ + 48// n + 32// + 16//n 3 - 40//n 2 - 232// n - 192// + 4//n 4 



64// 7/ 



• 48//n 2 + 388//n + 460// - 12//n 4 + 68// n 3 + 200//n 2 - 236// n - 558// 

2 ' ■■ 2 ■ ■ 16//// • 33//;r' • iW/ni 



-5//n 4 



4//n 3 



166/Zn 2 



15//n + 356// + 8/m 4 
- 110/i - 4n 4 + 6n 2 - 14n + 12)^C 2 (G) 
/[8(2/i -!)(// + n - l)(n - l) 3 (/i - 2)(n + 2)(n 2 - l)nT(i?)] 



(4.22) 
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Therefore summing all the contributions from each of the projections the full contribution from 
the final graph of fig. 2 is simply 

+ [32//n 2 + 32// n + 32 / u 7 - 96//V - 96^ 6 n - 160/x 6 + 28// n 4 + 56// n 3 - 20//n 2 

- 48// n + 324// - 148//n 4 - 296//n 3 + 394//n 2 + 542//n - 338// + 4//n 6 
+ 12//n 5 + 303//n 4 + 586//n 3 - 606// n 2 - 897//n + 188// - 8//n 6 - 24//n 5 

- 288//n 4 - 536//n 3 + 421//n 2 + 685//n - 50// + 5/m 6 + 15/m 5 + 125/m 4 
+ 225/m 3 - 144/m 2 - 254/m + Au - n 6 - 3n 5 - 19n 4 - 33n 3 + 20n 2 + 36n] 

x r(n + 1 - fj,)T(ji - l)ntfC 2 (G) 

/[8(2/x - 1)(m - 1) 2 (/" " 2)(n + 2)(n 2 - l)r(2 - n)T(ji + n + l)nT(R)} 

- (2/zn - n + l)r(/x - l)r(n)/if7i C 2 (G)/[2(2/i - 1)(M - 2 ) r (M + 

- [32//n 2 + 32// n + 32// - 144//n 2 - 144//n - 160// + 236// n 2 + 236// n + 316// 

- 4//n 3 - 172//n 2 - 160//n - 298// + 8/m 3 + 55 fin 2 + 31/m + 130/i - 4n 3 - 6n 2 

+ 6n - 20]/w7?C 2 (G)/[8(2/z - l)(/x - l) 3 (/z - 2)(n + 2)(n 2 - l)nT(i?)] (4.23) 

5 Results and discussion. 

We are now in a position to complete the calculation by adding the contributions from the four 
graphs of fig. 2 as well as including the gluon anomalous dimension in the Feynman gauge. Our 
final result is 

A+,i(a c ) = - [8//n 2 + 8//n + 8// + 2//n 4 + 4//n 3 - 22 //n 2 - 24// n 

- 28// - 6/m 4 - 12/m 3 + 14/m 2 + 20 fin + 32/t + 5n 4 
+ 10n 3 + n 2 - An - 12]r(n + 2 - /x)r(/x - l)fiC 2 {R)tf 
/[(u - 2) 2 (n + 2)(n + l)(n - l)r(2 - /i)r(/i + n)nT(R)} 

2u{y-l)S 1 {n)C 2 {G)r 1 ° 
(2»-l)(ji-2)T(R) 

- [32//n 2 + 32// n + 32// - 144//n 2 - 144//n - 160// - 4//n 4 

- 8//n 3 + 240//n 2 + 244//n + 316// + 16//n 4 + 32//n 3 

- 180//n 2 - 196//n - 306// - 20/m 4 - 40/m 3 + 59/m 2 
+ 79/m + 146// + 8n 4 + 16n 3 - 6n 2 - 14n - 28]//C 2 (G)r/? 
/[8(2/x - 1)0* - 1) 3 (^ - 2)(n + 2)(n + l)(n - l)nT(R)] 

+ [32//n 2 + 32// n + 32// + 8//n 4 + 16//n 3 - 120//n 2 - 128//n 

- 160// - 32// n 4 - 64//n 3 + 160// n 2 + 192// n + 316// + 48//n 4 
+ 96//n 3 - 78//n 2 - 126// n - 306// - 31/xn 4 - 62/m 3 + 31/m 

+ 14% + 7n 4 + 14n 3 + 7n 2 - 28]r(n + 2 - fi)T(fi - l)uC 2 {G)r,i 

/[8(2/i - 1)0* - 1) 2 (/* - 2)(n + 2)(n + l)(n - l)r(2 - n)F(ji + n)nT(J2)] 

(5.1) 

where we have set A + (a c ) = J2iZo ^+,i( a c)/Nf . One curious feature of this expression is that not 
all the T-function structures of the individual graphs survive. Specifically the terms involving 
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the factor T(fi — l)T(n)/T(fi + n) which appears in say ( |4.21 ) have cancelled in (5J). We can 



offer no explanation of this and it is not clear if it would be a feature of the analogous dimension 
in other (renormalizable) field theories. Moreover, the result is remarkably more compact than 



one would expect given the nature of (4,21) and (4.22) 



It remains to check that (|5.l| ) is in agreement with perturbation theory. Setting fj, = 2 — 
e and Taylor series expanding to several orders it is straightforward to produce the exact n- 
dependence for the two loop eigenvalue (|2.16| ). Moreover, expanding to three loops we obtain 
the following result for the combination of 3-loop entries of the mixing matrix, 

, &3ici _ 64(n 2 +n + 2) 2 (S 1 (n)) 2 C 2 (R) 
31 d n 3{n + 2)(n+l) 2 {n-l)n 2 T(R) 

64(10n 6 + 30n 5 + 109n 4 + 168n 3 + 155n 2 + 76n + 12)Si(n)C 2 (R) 
9(n + 2)(n + l) 3 (n - l)n 3 T(R) 

- 4[33n 10 + 165n 9 - 32n 8 - 1118n 7 - 5807n 6 - 12815n 5 - 16762n 4 - 13800n 3 

- 7112n 2 - 2112n - 288]C 2 ( J R)/[27(n + 2)(n + l) 4 (n - l)n 4 T(i?)] 
8(8n 6 + 24n 5 - 19n 4 - 78n 3 - 253n 2 - 210n - 96) Si (n)C 2 (G) 
27(n + 2)(n + l) 2 (n - l)n 2 T{R) 

- 2[87n 8 + 348n 7 + 848n 6 + 1326n 5 + 2609n 4 + 3414n 3 + 2632n 2 

+ 1088n + 192]C 2 (G)/[27(n + 2)(n + l) 3 (n - l)n 3 T(R)} (5.2) 

This needs to be compared with the first few non-zero moments of Q. To this end we have 
evaluated (|5.2| ) for even n, 2 < n < 24, and recorded the values as exact fractions in table 1. The 
coefficients with respect to each colour Casimir have been listed separately. The entries with n 
< 8 are in exact agreement with the explicit 3-loop MS expressions of || at this order in 1/Nf. 
Therefore we believe ( |5.1[ ) is correct. The provision of higher moments in table 1 is to be an aid 
to future checks on the explicit full analytic result as a function of n when it becomes available. 
Ordinarily in the perturbative calculation one has an additional check on the result emanating 
from general considerations. For the case when n = 2 the singlet operators correspond to the 
energy momentum tensor which is a conserved (non- anomalous) physical current. Therefore its 
anomalous dimension vanishes to all orders in perturbation theory, |2^] . Such a feature ought to 
be present in the large Nf results and this indeed was apparent in the dimension of the fermionic 
operator of [11]. However, in the gluonic case, as can be seen in table 1, the n = 2 entries are 
non-zero. This is not inconsistent with this principle. If one examines the combination of 
perturbative coefficients which vanishes at n = 2 one finds it is, at one loop, (a\ + b\) and (ci 
+ d\\Nf + d\2) with similar combinations at higher order. However the large Nf combinations 
are (a± — bic\/d\\) and (e?i 2 + b\Ci/dn) respectively for the fermionic and gluonic operators. 
Therefore clearly the former ought to be zero at n = 2 but not the latter. 



One feature of the form of the mixing matrix ( j2.4| ) is that the leading order series in Nf for 
7 99 (a) depends only on the Casimir C 2 (i?). This is clear from the one and two loop results of 
0, ^, |3|, |5| as well as the three loop results of Q for n < 8. Therefore if we make the assumption 
that this is true for all n at three loops then we can deduce the exact form of the coefficients 
c?3i in the C 2 (G) sector from (|5.2|). It is given by the last two terms of (5.2) which are each 



proportional to C 2 (G). Indeed if one examines the second column of table 1, one observes that 
the first four entries are the same as those which appear in ||, after noting that nf = 2T(R)Nf. 

As another motivation of this study is to provide a window into the structure of the operator 
dimension beyond what is currently available, it is a trivial exercise to expand (j5~J) for example 
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to four loops to deduce 

256(n 2 + n + 2) 2 [ff 3 (n) + 2(S 1 (n)) 3 ]C 2 (i?) 
27(n + 2)(n + l) 2 (ra - l)n 2 T(R) 

- 256[10n 6 + 30n 5 + 109n 4 + 168n 3 + 155n 2 + 76n + 12] 
x (S 1 (n)) 2 C 2 (R) / [27 (n + 2)(n + l) 3 (n - l)n 3 T(i?)] 

+ 256[37n 8 + 148n 7 + 697n 6 + 1573n 5 + 2087n 4 + 1725n 3 + 889n 2 
+ 264n + 36}S 1 (n)C 2 (R)/[81(n + 2)(n + l) 4 (n - l)n 4 T(i?)] 
256(n 2 + n + 2) 2 C(3)C 2 (i?) 
9(n + 2)(n + l) 2 (n - l)n 2 T(i?) 

- 16[77n 12 + 462n n + 1481n 10 + 3170n 9 + 12839n 8 + 37418n 7 
+ 75483n 6 + 103718n 5 + 98592n 4 + 64176n 3 + 27672n 2 
+ 7200n + 864]C 2 (.R)/[243(n + 2)(n + l) 5 (n - l)n 5 T{R)] 

32(9n 4 + 18n 3 + 79n 2 + 70n + 32)(Si(n)) 2 C 2 (G) 512Si(n)C(3)C 2 (G) 
+ 27(n + 2)(n + l) 2 (n - l)n 2 T{R) + 27T(i?) 

- 32[8n 8 + 32ra 7 + 95n 6 + 173n 5 + 781n 4 + 1311n 3 + 1176n 2 
+ 544n + 96]5i(n)C 2 (G)/[81(n + 2)(n + l) 3 (n - l)n 3 T(R)} 

_ 1024(n 2 + n + l)C(3)C 2 (G) 
27(n + 2)(n + l)(n-l)nT(i?) 

- 8[5n 10 + 25n 9 + 202n 8 + 658n 7 - 2639n 6 - 10115n 5 - 17896n 4 - 18216n 3 
- 11160n 2 - 3840n - 576]C 2 (G)/[243(n + 2)(n + l) 4 (n - l)n 4 T(i?)] 

(5.3) 

where £(n) is the Riemann (^-function. To obtain five and higher order combinations of pertur- 
bative coefficients at this order in 1/Nf one needs to first deduce the corresponding corrections to 
(|2.7| ) which are encoded in the 0(1/Nf) QCD /3-function, [H. Again if we make the assumption 
that 641 involves only C 2 (i?) then the terms proportional to C 2 (G) in ( |5.3D determine exactly 
the C 2 (Cr) part of d^. 

Finally we make some concluding remarks. First the determination of the expression ( |5.1| ) 
completes the leading order analysis in the large Nf expansion of the unpolarized anomalous 
dimensions for twist- 2 operators. To complete the full 0(1/Nf) analysis of deep inelastic scatter- 
ing one requires information on the corresponding coefficient functions. This has been achieved 
for the non-singlet case by various authors, |l5|, [TJ|]. As far as we are aware though, the un- 
polarized singlet case has yet to be examined. Second, it ought to be possible to extend our 
present work to determine the anomalous dimensions of the gluonic twist-2 singlet operators of 
polarized scattering. We would hope to return to this in a future article. 
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A Basic results. 



We give here several elementary integration rules for simple chain integrals. Although they have 
appeared before (see, for example, p3| ), we record the simplest cases here for completeness and 
also to fix our notation. First, we recall the chain result of for arbitrary a and j3, 



f 1 v(a,P,2n -a- (3) , 

h {y 2 ) a {{x-y?Y " (x 2 ) a +P-» { ' J 



where u(a,(3,"f) = ir^ ''a(a)a(ft)a(-f) and a(a) = T(/j, — a)/T(a) which is established by using 
Feynman parametrization. The more general moment type chain integral is given by 

f (Ay)P(A(x-y)) q . a . (Ax) p+q ,. . 

where i , m np{ a i^il) = ^ a m (a)a n (P)a p ('Y) , a n {a) = T(fi — a + n)/T{a) and m, n and p are 
integers. With flAj) one can, for example, build up two loop integrals which appear in the 
determination of the graph of fig. 3 where the exponent of the top left propagator is replaced 
by 0. 

Next we record the explicit expressions for the three loop recurrence relations we have used, 
where a and j3 are arbitrary and n is a strictly positive integer. Their derivation has been 
discussed in section 4. 

r ( R \ 2(2 l i + n){2^-a-f3 + n-2) ^ 

Cio 7 (a,/3,n) = — — — - — — Cio 7 (a + 1, /3, n + 2) 

(n + l)(n + 2) 

(2u — 2a + n — 3) „ . ^ _ , , 

(n + 1) 

+ [Cg£(a + 1, /?, n + 1) + C 007 (a, /?, n + 1) - + 1, /?, n + 1)] 

(2/*-a-/3 + n-2) 

X ; : 

(n + 1) 

- [O)o 7 (a + 1, A n + 2) + C$}(a + 1, P, n + 2) - C$(a + 1, /?, n + 2)] 
(2/i + n)(2/i-a-/? + ra-2) 



(n + l)(n + 2) 
[Cgg(a + 1, - 1, n + 1) - C™(a, p,n + l) 

+ Cg£(a + 1, /? - 1, n + 1) - C™(a, /3,n + l)]/(n + 1) (A.3) 



r « ^ 2(2^ + n+l)(2^-a-/? + n-2) ^ 

di 7 (a,/3,ra) = ? — — — — — Cn 7 (a + l,/3,n + 2 

(n + l)(n + 2) 

(2u — 2a + n — 3) _ . , „ „. 

- V , -Cn 7 a + l, /3, n + l 

(n + l) 

+ [2CgS} (a + l, /?, n + 2) - C o (7 -i) (a + l, /?, n + 2)] 
- a - + n - 2) 
X (n+l)(n + 2) 
+ [Cio 7 (a, /3, n + l) + Cg#(a + l, /3, n + l) - C7$}(a + l, A n + l)] 
(2^x - a - p + n - 2) 



(n + l) 



[Cio 7 (a + l, p, n + 2) + Cgi}(a + l, /?, n + 2) - C#°(a + l, /?, n + 2)] 
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(2/i + n + 2)(2/i -a-p + n-2) 
X (n + l)(n + 2) 

- 2[C 1 ° 1 2 7 (a + 1, /? - 1, n + 1) - Cg°(a, ft n + l)]/(n + 1) (A.4) 

r « ^ 2(2^ + n + l)(2^-a-/? + n-2) ^ 

C 207 (a,ftn) = , — — — — — C 2 o 7 (a + l,/3,n + 2) 

(n + lj(n + 2) 

(2ji-2a + n-3) _ . 

- , -tt -C 207 a + 1, A n + 1 

(n + lj 

2(2^-a-/? + n-2) 011 

+ ( n + l)(n + 2) W« + 1 >^ + 2 ) 

+ [C^Ja + l,ftn + 1) + Cio 7 (a, ft n + 1) - Cjg(a + 1, 0, n + 1)] 



x 



10 7 V" "T ? M) " T J. 7 T ^10 7 V", "Ti;- ^10 7 ^ 

(2fj, - a - + n - 2) 



(n + 1) 

- [Cio 7 (a + 1, ft n + 2) + CgJ}(a + 1, ft ra + 2) - Cj^(a + 1, ft n + 2)] 

(2/t + n + 2)(2/i - a- /? + n-2) 
X + l)(n + 2) 

- [Cgg (a + 1, - 1, n + 1) - C 2 12 °(a, ft n + 1) 

+ C °g(a + 1, /3 - 1, n + 1) - C 12 °(a, ft n + l)]/(n + 1) (A.5) 

We note that in their application to the integrals we require in this paper the first two terms on 
the right side which are the same function on the left but with a replaced by a + 1 will be zero 
by the dimension argument. 



B Explicit check for C 20 i(2/z - 4 + A, /i - 1 + A, 1). 

Throughout our discussion we have mentioned carrying out explicit checks on recurrence re- 
lations for low values of n. We give a detailed illustration of this here by summarizing the 
calculation of the result for C 2 Qi(2/i — 4 + A, fi — 1 + A, 1). Unlike the procedure used to con- 
struct recurrence relations we need to compute all of the invariant amplitudes which arise in the 
tensor decomposition of the corresponding integral without the null vector A^. At the end of this 
calculation contracting with the appropriate number of these vectors will give us a /^-dependent 
result which will agree with the general n-dependent expression evaluated at n = 1. Therefore 
if we label the indices of the y-propagator of fig. 7 as fi and v and that of the u-propagator as 
a, then the integral containing C 2 oi(2/i — 4 + A, \i — 1 + A, 1) decomposes into 

Arj^uXa + Blrj^Xu + rjuaXfj] + C ^-^-^ (B.l) 

The three contractions r]^ v x a , r}^ a x v and x^x v x a yield the respective combinations {2fiA + 2B 
+ C),(A+ (2fi + l)B + C) and (A + 2B + C). Using identities like xy = \[x 2 + y 2 - {x-y) 2 } 
we can rewrite the resulting integrals as a sum of scalar integrals. Again the majority of these 
are computable using the simple rules analogous to that of fig. 5. However, several basic results 



are necessary which it turns out have already been computed in the calculation of [14| and we 
merely quote these here 

r (0 l-LA (/i-l)V(l,l,2yu-2) 
Cooi(2^-4 + A,//- 1 + A,0) = 



r( M )A 

C 001 (2,-5 + A,,-l + A,0) = ,1,2^-2) (R2) 
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C 00 i(2/i-4 + A,^-2 + A,0) 



(ji - l) 2 (2/i i - 5/z + 4)i/ 2 (l, 1, 2// - 2) 



4I> + 1)A 
Therefore the above combinations become 

(2 M - 7)(p - 1)(M - 2)(/i - 3)i/ 2 (l, 1,2/* - 2) 



2[iA + 2B + C 
A + (2fi + l)B + C 
A + 2B + C 



6T(fi + 2)A 

(2/i 4 + ll// 3 - 121/i 2 + 374/x - 324) (p - l)i/ 2 (l, 1, 2/t - 2) 
48T(fi + 2) A 

(6/x 4 + 53// 3 - 333// 2 + 1408// - 1464) (/x - 1)^(1, 1, 2/i - 2) 



(B.3) 



96r(/i + 3)A 

Finally, solving for A, B and C the uncontracted integral is 
C 20 i(2//-4 + A,//-l + A,l) = 



7^x^(26/ - 261// 3 + 541// 2 - 576// + 120) 
+ (4/t 4 + 36// 3 - 143// 2 + 168// - 56) - 3) 

+ 2(2// 5 + 13/i 4 + 22/i 3 + 707// 2 - 1320// + 504) 



(//- 1)1/(1,1, 2/i -2) 



/[96(2/i-l)r(/i + 3)A] 



(B.4) 



whence the check on the recurrence relation we require is given by the coefficient of the x fM x u x a 
term. It is in agreement with the value of ([4.20 ) at n = 1. 
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n 


C 2 {R)/T{R) coefficient 


C 2 {G)/T(R) coefficient 


2 


_ 2716 

Q/1 Q 

z4o 


4232 

Q/1 Q 

z4o 


4 


3765671 
OU / oUU 


757861 
ou i ou 


6 


373918478 
1 zyoUo ( 


26390948 

OOOQ70 C 


8 


261337619387 

04UlUlozUUU 


420970849 

qoi A oonn 
oz!4oyUU 


10 


1994718278948 

QZuZOKJ < D4o / o 


2752314359 
zUo i OZ / yu 


12 


26908165162400891 

^ 7 n q o ? r; £; /i o q f; o n n 
DiVoZ ( ODD4ZODZUU 


2635361358193 
1 com oofin^/in 

ioyyiyzDyo4u 


14 


16449437567815379 


4616790551 

O.Z0U011UU 


16 


331641918853126821883 
f 028195 / 31 (539968000 


35631283224283 
24o83944o2o00 


18 


7462146988179117961999 
1578449269228939482375 


249465132233659 
16907213454300 


20 


8216273557478999257489 
1734325382898288180000 


1340260211455733 
89401072761000 


22 


228675683353467573965719 
48168621150965760797352 


648057715871521 
42614381628819 


24 


45507171465506238473037583 
9566587391661149850000000 


53358759817454947 
3463516806750000 



Table 1. Coefficients of \d%\ + bsici/du] as a function of moment. 
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Fig. 1. Operator insertion in gluon 2-point function. 



■e, 



Fig. 2. Leading order diagrams for A + (a c ). 




Fig. 3. Definition of G pqn (a). 
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Fig. 4. Definition of H pn (a,{3). 




Fig. 5. Example of elementary integral rule. 




Fig. 6. Intermediate integral in calculation of H pn {a,j3). 




Fig. 7. Definition of C pqi {a, j3,n). 
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